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Abstract 

In this paper, the stabihzed finite element method based on lo- 
cal projection is applied to discretize the Stokes eigenvalue problems 
and the corresponding convergence analysis is given. Furthermore, 
we also use a method to improve the convergence rate for the eigen- 
pair approximations of the Stokes eigenvalue problem. It is based on 
a postprocessing strategy that contains solving an additional Stokes 
source problem on an augmented finite element space which can be 
constructed either by refining the mesh or by using the same mesh 
but increasing the order of mixed finite element space. Numerical 
examples are given to confirm the theoretical analysis. 
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1 Introduction 

In this paper, we are concerned with the Stokes eigenvalue problems. The 
study of Stokes eigenmodes is required when the dynamics behaviors gov- 
erned by the Navier-Stokes equations result from the way this nonlinear dy- 
namics is controlled by diffusion. For the other reasons to study the Stokes 
eigenmodes, please read the papers [6| [22]. 
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The Stokes eigenvalue problem reads as follows: 
Find (u,p, A) such that 



Au + Vp 
Vu 



Au in 

in 

on dn 
1, 



(1.1) 



u 



In ^'d^ 



where C is a bounded domain with Lipschitz boundary dQ and A, V, 
V- denote the Laplacian, gradient and divergence operators, respectively. 

There are several works for the eigenvalue problems and their numerical 
methods such as Babuska and Osborn [21 [3], [29], Mercier, Osborn, Rappaz 
and Raviart [27], etc. Osborn [29], Mercier, Osborn, Rappaz and Raviart [ 27 \ 
give an abstract analysis for the eigenpair approximations by mixed/hybrid 
finite element methods based on the general theory of compact operators 
([13j). In [2T] and [21], a posteriori error estimates and the corresponding 
adaptive finite element methods are given for the Stokes eigenvalue problems. 

The first aim in this paper is to use the focal projection stabilization 
method to discretize the Stokes eigenvalue problems. The local projection 
stabilization (LPS) method has been proposed for the Stokes problem in [7]. 
The extension to the transport problem was given in [8] . The analysis of the 
local projection method applied to equal-order interpolation discretization 
can be found [25] for Oseen problem and in for convection-diffusion 
problem. 

The stabilization term of the local projection method is based on a pro- 
jection TT/i : \4 — 7- Vfi of the finite element space Vh which approximates the 
solution into a discontinuous space V^. The standard Galerkin discretiza- 
tion is stabilized by adding a term which gives control over the fiuctuation 
id — 7ih of the gradient of the solution. Here, the LPS method is based on the 
approximation space Vh and the projection space are defined on the same 
mesh. In this case, the approximation space Vh is enriched to satisfy the lo- 
cal inf-sup condition guaranteeing the existence of an interpolation with an 
orthogonal property compared to standard finite element spaces. For more 
details, please read the book [31|- 

Recently, many effective postprocessing methods that improve the con- 
vergence rate for the approximations of the eigenvalue problems by the finite 
element methods have been proposed and analyzed ([HISniES])- Xu and Zhou 
[35] have given a two-grid discretization technique to improve the convergence 
rate of the second order elliptic eigenvalue problems and integral eigenvalue 
problems. Racheva and Andreev [30], Andreev, Lazarov and Racheva [1] 
have proposed a postprocessing method that improve the convergence rate 
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for the numerical approximations of 2m-order selfadjoint eigenvalue prob- 
lems especially biharmonic eigenvalue problems. In [13], a similar method 
has been given for the Stokes eigenvalue problem by mixed finite element 
methods. The second aim of this paper is to propose and analyze a postpro- 
cessing algorithm which can improve the convergence rate of the eigenpair 
approximations for the Stokes eigenvalue problem by the LPS method. 

The postprocessing procedure can be described as follows: (1) solve the 
Stokes eigenvalue problem in the original finite element space; (2) solve an 
additional Stokes source problem in an augmented space using the previous 
obtained eigenvalue multiplying the corresponding eigenfunction as the load 
vector. This method can improve the convergence rate of the eigenpair ap- 
proximations with relative inexpensive computation because we replace the 
solution of the eigenvalue problem by an additional source problem on a finer 
mesh or in a higher order finite element space. 

An outhne of the paper goes as follows. In Section 2, we introduce the 
application of LPS method for Stokes eigenvalue problem. The corresponding 
error estimate is given in section 3. Section 4 is devoted to deriving the 
postprocessing technique and analyze its efficiency. In Section 5, we propose a 
practical computational algorithm to implement the postprocessing method. 
In Section 6, we give two numerical results to confirm the theoretical analysis. 
Some concluding remarks are given in the last section. 



2 Discretizations of the Stokes eigenvalue prob- 



In this paper, we use the standard notations ([HI [T^l [IS]) for the Sobolev 
spaces H'^{Q) (standard interpolation spaces for real number m) and their 
associated inner products (■, norms || ■ \\m and seminorms | ■ \m for m > 0. 
The Sobolev space H^{Q) coincides with L^(f2), in which case the norm and 
inner product are denoted by || ■ || and respectively. In addition, denoted 
by Lq{Q) the subspace of L^(f2) that consists of functions on L'^{Q) having 
mean value zero. We also use the vector valued functions {H"^{Q)y just as 
[12] and [20]. 

The corresponding weak form of flLip is: 

Find {u, p, X) e V X Q X R such that r(u, u) = 1 and 



lem 



a(u, v) - b{-v,p) 
b{u,q) 



Ar(u,v) Vv G V, 
Vg G Q, 



(2.1) 
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where V = {H^{n))\ Q = Ll{n) and 

a(u,v) = / VuVvrffi, 
Jn 



uv(if2. 

From we know eigenvalue problem (12.1 p has an eigenvalue sequence {Aj}: 
< Ai < A2 < • ■ ■ < Afc < ■ ■ ■ , lim Afc = 00, 

fc— >-oo 

and the associated eigenfunctions 

(Ui,pi), (U2,P2), ■ ■ ■ , (Ufc,Pfc), ■ ■ ■ , 

where r(uj,Uj) = 5ij. 

For the aim of analysis, we define the bilinear form as 

A((u,p);(v,g)) = a(u,v)-6(v,p) + 6(u,g). (2.2) 

For simplicity, we only consider the simple eigenvalues in this paper. We 
know that a(-, ■), 6(-, ■) and s(-, ■) have the following properties ([20]): 

a(u,v) < ||u||i||v||i, (2.3) 

a(u,u)>C||u||2, (2.4) 

r(u,v)<C||u||o||v||o, (2.5) 

r(u,u)>C||u||2, (2.6) 

sup ^-p^>C\\q\\o. (2.7) 
Ot^vEV ||v||i 

||u||i + IIpIIo < C sup —-- — —-- — , (2.8) 

07^(v,g)eVxQ ||V||l + ll^llo 

where C > 0. In this paper, C denotes constant independent of the mesh size 
h and sometimes depends on the eigenvalue A and may be different values at 
its different occurrence. 

For the eigenvalue, there exists the following Rayleigh quotient expression 

a(u, u) 

A= ) ' ! . (2.9) 
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2.1 Local projection stabilization 

In this section, we consider equal order interpolations stabilized by the local 
projection method in its one-level variant as developed in [191 ES]. For the 
two- level approach we refer to [71 HOl |28]. Let Vh denote a scalar finite 
element space of continuous, piecewise polynomials over 7/,. The spaces 
for approximating velocity and pressure are given by \h '■= l~l ^ 
Qh '■= Vh n Q. The discrete problem of our stabilized method is: 

Find {uh,ph) eVhX Qh such that 

f a{uh,v) - b{v,ph) = Xhr{uh,v) Vv G V/,, 
\ b{uh,q) + Sh{ph,qh) = yqeQh, 

where the stabilization term with user-chosen parameters ax is given by 

Sh{p,q) = aK{i^h'^P,Hh'^q)K- (2.11) 

KeTh 

Here, the fiuctuation operator Kh : L^(f2) — )• L^(f^) acting componentwise is 
defined as follows. Let Ps{K) denote the set of all polynomials of degree less 
than or equal to s and let Dh{K) be a finite dimensional space on the cell K G 
Th with Ps{K) C Dh{K). We extend the definition by allowing P-i{K) = 
Dh{K) = {0}. We introduce the associated global space of discontinuous 
finite elements 

Dh := DhiK) 

KeTh 

and the local L^(K)-projection ttk '■ L'^{K) — > Dh{K) generating the global 
projection tt/j : L^(f2) — )■ Dh by 

{'Khw)\j^:= 'Kk{w\k) WK eTh, e L^{Q). 

The fiuctuation operator Kh : L'^{Q) — t- L'^{Q) used in (12. lip is given by 
Kh '■= id — TCfi where id : L^(fi) — > L'^{Q) is the identity on L'^{Q). 

In order to study the convergence properties of this method for Stokes 
eigenvalue problem, we introduce the bilinear form 

Ah{{u,p); (v, q)) = (Vu, Vv) - {p, div v) + {q, div u) + Sh{p, q). (2.12) 

and the mesh- dependent norm 

|||(v,^/)||U:= (|v|?+ 11^11^+ 5^ «x||/^.Vg||^,A') • (2.13) 
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The existence and uniqueness of discrete solutions of Stokes problem have 
been studied in [251 HI] for different pairs {Vh,Dh) of approximation and 
projection spaces, respectively. Based on these results, the existence and 
uniqueness of eigenvalue problem fl2.10p can be given similarly. 

The stability and convergence properties of the LPS method fl2.10p need 
the following assumptions ( [2^1 EI])- 

Assumption Al: There is an interpolation operator : H'^{Q) — t- Vh such 
that 

11^^ - 'ihv\\o,K + hxlv - ihv\i,K < Ch[\\v\\i^i^(_K) (2.14) 

for all K E Th, v G H\u{K)) and 1 < / < /c + 1, where oj{K) denotes a 
certain local neighborhood of K which appears in the definition of these in- 
terpolation operators for non-smooth functions; see [T71 [33] for more details. 
Assumption A2: The fluctuation operator satisfy the following approx- 
imation property 

\MWk < Ch'j,\q\i^K yK Erh,yqeH\K),0<l<k. (2.15) 



Assumption A3: There exists a constant /3i > such that for all /i > 

inf sup ,, ii^'^^'if^i > /3i > (2.16) 

qh&T^h{K) .i,^^Vh{K) \\Vh\\o,K \\<lh\\o,K 

is satisfied where Vh{K) = {vh\K '■ Vh E Vh, Vh = in Q\K}. 

The assumption Al and A3 guarantee the existence of an interpolant with 
the usual interpolation properties (12.141) and the orthogonality 

iv-JhV,qh) = yq^E V^^^v E H\9), (2.17) 

whereas A2 is needed to bound the consistency error pT]. For example, in 
the one-level LPS assumption Al and A2 are satisfied if we choose (V/i, Dh) = 
(Pfci-Pfci^i) continuous and discontinuous, piecewise polynomials of degree r 
and r — 1, respectively. In order to guarantee A3, Vh need to be enriched by 
suitable bubble functions. For more details about LPS method, please read 
the papers [25l [26] and the book [3T] . 

Lemma 2.1. (flMI) Let the assumption Al, A3, and ax ~ h'j^ be fulfilled. 
Then, there is a positive constant j3A independent of h such that 

■ r ^h((vh,gfe);(wfe,rfe)) 

mf sup --T- --r^ — — > Pa>0 (2.18) 

holds. 

Based on Lemma 12.11 the discrete Stokes eigenvalue problem (12.101) is 
consistent with the continuous problem (12. ip ([19j). 
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3 Convergence analysis 



In thois section, we give the convergence analysis for the eigenpair approxi- 
mation {uh,Ph,)^h) in fl2.10p . 

We know that the convergence rate of the eigenpair approximations by 
the finite element methods depends on the regularities of the exact eigen- 
functions. The exact eigenfunctions of the Stokes problem only belong to 
the space {H^{Q)y x H^[Q) on general domains. But for the domains with 
smooth boundary, the exact eigenfunctions have additional regularities. In 
this case we need to use isoparametric mixed finite element methods to fit 
the domain more exactly ([TT] and [IS])- The goals of this paper are to use 
LPS method to solve the Stokes eigenvalue problem, and propose and an- 
alyze a postprocessing method which can improve the convergence rate for 
both eigenvalue and eigenfunction approximations. The assumption that Q 
is a convex polygonal domain can make the expression of the main idea of 
this paper more directly. But, we need to notice that this assumption limits 
the regularity of the exact eigenfunctions and makes the analysis of the con- 
vergence rates much more complicated. It is well known (|H|5l[T8]) that for 
a given f G {H'^{Q))'^ the solution (u,p) of the corresponding Stokes problem 

a(u, v) - 6(v,p) = r(f,v) Vv G V, , . 

b{u,q) = yqeQ ^ ' 

has the following regularity ( [H El El [IE] ) 

||u||2+^ + |b||i+-, < C||f||^ Wie{H\n)f, (3.2) 
where < 7 < 1 is a parameter that depends on the largest interior angle of 

dn (i). 

From fl2.10p . we can know the following Rayleigh quotient for A/, holds 

a{uh,Uh) 

A/i - r [6.6 

Holds. 

It is also known from [3] the Stokes eigenvalue problem (12.101) has eigen- 
values 

< (Ai)^ < iX2)h <■■■< i\k)h <■■■< {\N)h, 
and the corresponding eigenfunctions 

((Ul)/i, (Pl)ft), (Mh, iP2)h), • ■ ■ , ((Ufc)/i, iPk)h), ■ ■ ■ , ((uat)/!, {pN)h), 

where r((uj)/i, (uj)/i) = 6ij,l < i,j < N, N denotes the dimension of the 
finite element space Yh x Qh ■ 
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Let us define the compact operator T : (L^(fi))^ — )■ {H^{fl)y and the 
operator K : (^^(fi))^ ^ L2(fi) by 

A((Tf,Kf),(v,g)) = r(f,v), V(v, G x ^^(fi). (3.4) 

Hence the eigenvalue problem fl2.ip can be written as 

ATu = u. (3.5) 

Let M(Aj) denote the eigenspace corresponding to the eigenvalue Aj which is 
defined by 

M{Xi) = ^^{w,tlj) e {H^{Q)y X Ll{Q) : (w. tb) is an eigenfunction of fl2J|) 
corresponding to Aj and r(w, w) = l|- 

Similarly, we also introduce the discrete operator Th : (L^(fi))^ — t- Yh 
and the operator Kh : (L^(i7))^ — )■ Qh by 

Ah{{Thf,Khf),{v,q)) = r(f,v), V(v, q) G V;, x Q;,. (3.6) 

Hence the operator form of the discrete eigenvalue problem fl2.10p is 

XhThUh = Uh. (3.7) 

In [T9], the convergence result of LPS method for Stokes problems has 
been given. Combining abstract spectral approximation results from [3], we 
can give the convergence results for the Stokes eigenvalue problem by LPS 
method. The eigenvalue approximation and the corresponding eigenfunc- 
tion approximation {uh,ph) have the following error bounds ([21 [2H1 [ISl [13 
EZIEO]): 

||u-U;,||l < C||(T-T/,)|m(a)||i. (3.8) 

In order to do the analysis for the postprocessing in the following sections, 
we also need the convergence result for the eigenfunction approximation Uh 
in iJ~^-norm. For this aim, based on the result in [2], we first need to use 
the duality argument to get iJ~^-norm error estimate for the finite element 
projection and the process is similar to the one in the paper for the 
L^-norm error estimate. 

The finite element projection {RhU, Rhp) denotes the finite element solu- 
tion of the following Stokes problem: 

Find {Rh{u,p),Gh{u,p)) G V/i x Qh such that 

Ah{{Rh{n, p), Gh{u, p)); (v, q)) = A{{u, p); (v, q)) V(v, g) G V;, X Q,. (3.9) 
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(3.10) 



||u-u,||i < C\\{T-Rf,{T,K)\Mix)\\i. (3.11) 

Prom the definition, we have the orthogonal relation 

Ah{{u- Rh{u,p),p-Gh{u,p));{v,q)) = Sh{p,q) 

V(v,g)eV^xgft.(3.12) 



|||(vft - Rh{n,p),qh - ^^(u,^))!!^ 
^ J_ Afe((vfe - Rh{u,p),qh - Gh{u,p)); (wfe,^;^)) 

~ Pa o^{y,h,iPh)(^VhxQh |||(w/j,V'/i)||U 

< J_ Afe((vfe-u,gfe-p);(wfe,V;fe)) 

~ 0^(w^,Vfe)GVhxQ^ |||(Wfe,V'/i)||U 

_^J_ A((vi - Ri,,{u.p),p - Gu{u,p))-, (w/,,, ( vj) 

Pa o^i^^hMeVhxQh \\\{^h,'4^h)\\\A 

< -3-\\\{^-^h,P- qh)\\\A + ^ sup 



/^A 07^(wfe,Vft)eVaxQ^ |||(wft, V^/i)||U 

< ^\\\{n-v,,p- qn)\\\A + ^Sl/\P,P) (3.13) 
Pa pa 



|||(u - i?,,(u,p),p - G/,(u,p))||U 

< lll(u - Vh,p - g/,)||U + |||(vft - Rh{u,p),qh - ^^(u,^))!!^ 

< C\\\{u-VH,p-qh)\\\A + CS'J'{p,p). (3.14) 

The arbitrariness of (v/^, qh) leads to the following 

\\\{n-Rh{n,p),p-Gh{u,pm\A 
< G{ inf |||(u-v,,p-g,)IIU + ^f (P,P)). (3.15) 

(vft,gh)6VhxQft 

||u-u^||i < C5ft(A), (3.16) 

where 5/j(A) is defined by 
Sh{X):^ sup ( inf \\\{w-Vh,^-qhmA + S'J'{^,i;))i3.17) 

(w,V')gM(A) (vft,9h)eVhXQh 
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We choose g G {H^{n)) 



^ such that II g 111 



1 and 



u - Rh{u,p)\\ 



r(u - Rh{u,p),g). 



Then we define a duahty problem corresponding to g: 
Find {ug,pg) G V x Q such that 



Combination of (13.121) and (13.181) derives the following estimate 
r{u- Rh{u,p),g) = a{u- Rh{u,p),Ug) -b{u- Rh{u,p),pg) 



= a(u - Rh{u,p),Ug - \h) - b{u - Rh{u,p),pg - qt) 
-b{Ug - \h,p-Gh{\y,p)) - ShiGhin,p),qh) 

= a(u - Rh{u,p),Ug - v/J - 6(u - Rh{u,p),pg - qh) 
-b{ug - ^h,p- Gh{u,p)) - Ship - Gh{u,p),pg - Qh) 

+ Ship- Ghiu,p),Pg) + Sh{p,Pg - Qh) - Ship,Pg). 



Choosing iyhAh) ^ x Qh as an interpolant of (ug,pg), we obtain 



|r(u-i?;,(u,p),g)| 

< C(||u - i?,(u,p)||i + lb - G',,(u,p)||o)(||u, - v;,||i + Ib^ - g,,||o) 
+ \Sh{Gh{}i,p),qh)\ 

< C(|||(u - Rh{m,p).P- G,,(u,p))||U + s];\p,p)) 

{\\\{Ug-Wh,Pg~ qh)\\\A + S]l^{Pg,Pg)). 



In particular, when Vt is smooth, we have the regularity estimate 



«(v,Ug) - b{w,pg) 



r(v, g) Vv G V, 
Vg G Q. 



(3.18) 



U3II3 + Ibglh < C||g||l, 



U - i?;,u|| 



< Cr]h5h 



where 



rjh= sup inf {\\\{Tg-^r,Kg-q)\\U + Sl/\Kg,Kg)). (3.19) 

||g||l = l {v,q)&VhXQh 
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4 One correction step 



In this section, we present a type of correction step to improve the accuracy 
of the current eigenvalue and eigenfunction approximations. This correction 
method contains solving an auxiliary source problem in the finer finite ele- 
ment space and an eigenvalue problem on the coarsest finite element space. 
For simplicity of notation, we set (A, u) = (Aj, Ui) {i = 1,2, k, ■■ ■) and 
{Xh,Uh) = {Xi,h,Ui^h) {i = 1,2,- ■■ , Nh) to denote an eigenpair of problem 
(??) and (??), respectively. 

To derive our method, we need first to introduce the error expansions of 
the eigenvalues by the Rayleigh quotient formula. It is well known that there 
have been the Rayleigh quotient error expansions for the eigenvalues of the 
second order elliptic problems (|23]). 

Theorem 4.1. Assume (u,p, A) is the true solution of the Stokes eigenvalue 
problem i\2l\} . w e {H^{Q))^ and i/j G Ll{Q) satisfy 



Let us define 



b{w,i:) + SH{iJ,ij) = 0. (4.1) 



A =4^. (4.2) 
r(w, w) 



Then, we have 

- a(w — u, w — u) + 2b{w — u, p — -0) — Ar(w — u, w — u) 

A - A = 

r(w, wj 

25h(^,^) 



(4.3) 



r[w, w) 

If the condition i| ) is changed to be 

b{w,^) = 0, (4.4) 

the expansion for A — A should be 

- a(w — u, w — u) + 2b(w — u, j9 — -0) — Ar(w — u, w — u) 

A-A = ^ . 4.5 

r(w, wj 

Proof From ([2A]), fl2A0|l . flOjl . (gj]), (g2]) and direct computation, we have 

- a(w, w) — Ar(w, w) 

A - A = 

r w, w 



11 



a(w — u, w — u) + 2a(w, u) — a(u, u) — Ar(w, w) 
r(w, w) 

a(w — u, w — u) + 2Ar(w, u) + 2fo(w, p) — Ar(u, u) — Ar(w, w) 

r(w, w) 

a(w — u, w — u) — Ar(w — u, w — u) + 26(w, p) 
r(w, w) 

a(w — u, w — u) — Ar(w — u, w — u) + 26(w, p — i') — 2Sh{'ip, "0) 

r(w, w) 

a(w — u, w — u) — Ar(w — u, w — u) + 26(w ~ u,p — ip) — 2Sh{'ip, "0) 



This is the desired result (14.31) and the expansion of (14.51) can be prooved 
similarly. □ 

Assume we have obtained an eigenpair approximation (A^j, u/t^,p/ij G 
71 X Yfii X Qh- Now we introduce a type of correction step to improve the 
accuracy of the current eigenpair approximation (A/^^, u/ij,p/ij. Let V/^j x 
Qh2 C (i^Q (fi))^ X Ll{fl) be a finer finite element space such that V/^ x Q^^ C 
Vftj X (5^2- Based on this finer finite element space, we define the following 
correction step. 

Algorithm 4.1. One Correction Step 

1. Define the following auxiliary source problem: 
Find ^ such that 

a{uh2,Vh2) = \h^b{uh^,Vh2), Mvh^eVh^. (4.6) 

Solve this equation to obtain a new eigenfunction approximation G 

Vh2- 

2. Define a new finite element space VH,h2 = Vh + spanlw/t^} and solve 
the following eigenvalue problem: 

Find {Xh2,Uh2) G 7^ x VH,h2 such that b{uh2,Uh2) = 1 and 

0'{Uh2,VHM) = ^h2b{Uh2,VHM)^ ^VH,h2 ^VH,h2- (4-7) 

Define (A^ajW^J = Correction{VH, Xhi,Uh^,Vh2)- 

Theorem 4.2. Assume the current eigenpair approximation {Xh^, Uh^) eTZx 
V/ii has the following error estimates 

\\u-Uh,\\a < £fti(A), (4.8) 
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\\u-Uh,\\-a < r]a{H)\\u - Uh^Wa, (4.9) 

lA-A.J < 4(A). (4.10) 

Then after one correction step, the resultant approximation {Xh^iUh^) eTZx 
V/i2 has the following error estimates 

\\u-UhJa < eh2(A), (4.11) 

\\u-Uh.2\\-a < r]a{H)\\u - Uh^Wa, (4.12) 

|A-A,J < 4(A), (4.13) 
where en^iX) := 7]a{H)eh^{X) + el^{X) + ^^^(A). 

Proof From problems (??), (??) and (jM]), and and fHTTU]) . the 

following estimate holds 

\\uh2 - Ph2u\\l ^ a{uh2 - Ph2U,Uh2 - PhiU) = b{Xh,Uh, - Xu,Uh2 - PhiU) 

< ||A/jiM/ii - AM||_a||M/i2 - -P?»2w||a 

< (|Ah,, - A|||Mh,J|_a + X\\Uh^ - u\\^a)\\Uh2 " ^ha^lU 
^ (4i(A) +r/a(^)e/.i(A))||Sft2 -Pft2«||a. 

Then we have 

\\uh2-Ph2u\\a < el{X) + ria{H)ehA^). (4.14) 
Combining fl4.14|) and the error estimate of finite element projection 

\\u-Ph2u\\a < <5/,2(A), 

we have 

\\uh2-u\\a < 4(A)+r7„(if)£^,(A) + 5;,,(A). (4.15) 

Now we come to estimate the eigenpair solution (A/ijjtt/jJ of problem fl4.7l) . 
Based on the error estimate theory of eigenvalue problem by finite element 
method ([?,|3]), the following estimates hold 

Wu-Uh^Wa < sup inf \\W - V\\a < \\U - Uh^Wa, (4.16) 

and 

\\u-Uh2\\-a < rja{H)\\u - Uh2\\a, (4.17) 

where 

^a{H) = sup inf \\Tf -v\l<ri,{H). (4.18) 

From fHl5|) . firTHjl . fHTTD and fliTTSjl . we can obtain fl4llD and fHA2|) . The 
estimate fl4.13p can be derived by Theorem ?? and (14.111) . □ 
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If the eigenpair approximation {uh,ph, Xh) of the Stokes eigenvalue prob- 
lem f l2.ip has been obtained, we define the following Stokes source problem: 
Find (u,p) G V X Q such that 

a(u,v) - 6(v,p) = Xhr{uh,v) Vv G V, . - 

6(u,g) = VgGQ. ^ ' 

We also define the following Rayleigh quotient formula for the solution {u,p) 

~ a(u, u) 

A = 4.20 
r(u, u) 

For the eigenpair {u,p, A), we can give the following error estimate. 

Theorem 4.3. Assume (u,p, A) is the true solution of the Stokes eigenvalue 
problem i\2. (u/j,]?^, A/j) is the corresponding finite element solution of the 
discrete Stokes eigenvalue problem Ii2.10\) . (u,]5) is the true solution of prob- 
lem Iji4-^^ '^iT'd A is defined by ^^^W^. Then we have the following estimates 

||u-ii||i + IIp-pIIo < C(||u-u/,||_i + |A- A,,|), (4.21) 
|A - A| < C(||u - VLh\ti + |A - \h?). (4.22) 

Proof. First from Stokes eigenvalue problem fl2.ip and Stokes problem fl4.19p , 
we have 

a(u - u, v) + 6(v, p -p) + 6(u - u, q) 
=r{\h\i.h - Au, v) 
=\hr{\ih - u, v) + (A/, - A)r(u, v) 

<C(||u;,-u||_i + |A/,-A|)||v||i. (4.23) 
Then, from f l2.8p . we have 
II ~ II , ii~ II / a(u-u,v) + 6(v,p-j9) + 6(u-u,g) 

||U — U||i + Hp — pIIo < sup T — T 7—r 

07^(v,g)GVxQ ||v||l + \\q\\o 

<C(||u,-u||_i + |A,-A|). (4.24) 



From fl4.24p and the Rayleigh quotient expansion (14. 5p . we obtain ([IS]) 

A-A < C(||u - u||^ + ||u - u||i||p - pIIo) 

< C(||u^-u||_i + |A;,-A|)2 

< C(||u,-uf_^ + |A-A,|2). 

So the proof is complete. □ 
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Based on the result of the convergence rate for the eigenpair approxima- 
tion, we can obtain the error estimates: 

For the smooth domain, from (??)-(??) and (??)-(??) 

||u- u||i + ||p-jo||o < C/i^ forA; = l, (4.25) 

\~X-X\<Ch^ {oTk = l, (4.26) 

\\u-u\\i + \\p-p\\o<Ch''+'^ for A; > 2, (4.27) 

|A - A| < C7/i2fc+^ for k>2. (4.28) 

For the convex polygonal domain, from (??)-(??) and (??), we have 

\\u-u\\i + \\p-p\\o<Ch^', (4.29) 
\~X-X\<Ch^'. (4.30) 

This means that (u,p, A) is a better approximation than {uh,Ph,Xh) of the 
true solution {u,p, A) of the Stokes eigenvalue problem (12. ip . 

5 Multi-level correction scheme 

In this section, we introduce a type of multi-level correction scheme based on 
the One Correction Step defined in Algorithm 14.11 This type of correction 
method can improve the convergence order after each correction step which 
is different from the two-grid method in |35] . 

Algorithm 5.1. Multi-level Correction Scheme 

1. Construct a coarse finite element space Vh and solve the following 
eigenvalue problem: 

Find {Xh, uh) xVh such that hiuu, uh) = 1 and 

a{uH,VH) = XHb{uH,VH), ^Vh&Vh- (5.1) 

2. Set hi = H and construct a series of finer finite element spaces Vh^, ■ ■ ■ ,Vh, 
such that r]aiH) > ShA^) > 6h,iX) > ■ ■ ■ > 5hM- 

3. Do k = 0,l,--- ,n-2 

Obtain a new eigenpair approximation (A/^^.^^, ■U/^^.^J E 71 x Vh,,^^ by a 
correction step 

(^/^fc+i'^/ifc+i) = Correction(VH, Xh^,Uh^,Vh^^J. (5.2) 

end Do 
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4- Solve the following source problem: 
Find Uh„ G Vh„ such that 

a{uh„,VhJ = \h„_Xvhr,^i,Vh„), ^Vh„eVh„. (5.3) 

Then compute the Rayleigh quotient of u^^ 

Finally, we obtain an eigenpair approximation {\hn,UhJ) ElZx Vh„- 

Theorem 5.1. After implementing Algorithm \5.1\ the resultant eigenpair 
approximation {Xh^,Ufi^) has the following error estimate 

\\uh,,-u\\a < £hM, (5.5) 
|A.„-A| < eliX), (5.6) 

n 

where £/,„(A) = ^ ^a(^/)"-'5h,(A). 

k=l 

Proof From r]a{H) > Sh,{X) > ^/.^(A) > ■ ■ ■ > and Theorem S^l we 

have 

^h.+.(A) < r^a(ff)£,,(A) + 5,,^,(A), forl<A;<n-2. (5.7) 
Then by recursive relation, we can obtain 

£^„-.(A) < VaiH)ef,,U>^) + 6,^_AX) 

< Va{H)'e,„_,{X) + r^a{H)S,„_,{X) + 5,_,(A) 

n— 1 

< Y.'^aiHr-'^-%,{X). (5.8) 

k=l 

Based on the proof in Theorem 14.21 and (15.81) . the final eigenfunction approx- 
imation Uh^ has the error estimate 

\\uh„-u\\a < £l^_^iX) + riaiH)eh,,_,{X) + 5h,XX) 

n 

< Y.'^a{Hr-%SX). (5.9) 

fc=l 

This is the estimate (15.51) . From Theorem ?? and (15. 9p . we can obtain the 
estimate (15. 6p . □ 
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6 Postprocessing algorithm 



Theorem 14.31 has only theoretical value and cannot be used in practice since 
the exact solution of the Stokes source problem fl4.19p is always not known. 
In order to make it useful, we need to get a sufficient accurate approximation 
of the Stokes source problem. Here we discuss two possible ways how to 
obtain the approximation of the Stokes source problem f l4.19p . The first way 
is the so-called "two-grid method" of Xu and Zhou introduced and studied 
in [35] for second order differential equations and integral equations. The 
second way proposed and studied by Andreev and Racheva in [30] uses the 
same mesh but higher order finite element space. 

The first way uses a finer mesh (with mesh size h'^) to get an approxima- 
tion of A with an error 0{h^^) or 0{h'^'^) for k < 2. The advantage of this 
approach is that it uses the same finite element spaces and does not require 
higher regularity of the exact eigenf unctions. The second way is based on 
the same finite element mesh Th but using one order higher finite element 
space. Also, to get an improvement for the approximation of Xh to the error 
0{h^) or Oih?^'^"') from O^h"^), we need to investigate the regularity of the 
Stokes eigenvalue problem. 

We can treat both ways in the same abstract manner. Namely, let us 
introduce the enriched finite element space x such that Yh ^ Qh 

^ Qh {Hl{yt)Y X Lliyt) and consider the following discrete Stokes 
problem: 

Find (u/i,pft,) G V/j x Qh such that 



Here, we suppose that the approximation (u^, ph) G V/^ x Q^ has the following 
error estimate: 

For a smooth domain 



So, we need define the following Rayleigh quotient for (u/i,pft) 

~ a{uh,^h) 

Xh = -7-~ — 6.4 

From the analysis above, we can obtain the following error estimate for 
the new eigenpair approximation (ufe,p/i,A/i) G V/i x x R. 




(6.1) 



||U - Uftlli + lb -Phllo < C'^''^^(||u||a:+2 + lb||fe+l) 

and for a convex polygonal domain 



(6.2) 



u - u/,||i + lb -Phllo < CY?Wyi\s+2 + Iblls+i)- 



(6.3) 
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Theorem 6.1. Assume Xh is defined by Ili6.4\ ), {vLh,Ph) is the solution of Ii6. 1\) 
and (u,p, A) is the true eigenpair of the Stokes eigenvalue problem 112.1]) . 
Then we have 

\Xh - A| < C(||u - Uh\\-i + |A - Xh\ + ||u - Uh\\i + \\p-Ph\\of 

+ CShiPh,Ph), (6.5) 

\\uh - u||i + \\ph -p\\o < C(||u - Uh\\-i + |A - Xh\ + ||u - Uhlli 

+ \\p-ph\\o)- (6.6) 

Proof. First from fl4.21l) and the triangle inequahty, we can obtain (16.61) . 
Using b{uh,ph) + Sh{ph,Ph) = and (14.31) . the foUowing error estimate holds 

|Aft-A| < C{\\uh-u\\l + \\ph-p\\l) + Sh{ph,Ph) 

< C(||u - Uh\\-i + |A - Xh\ + ||u - UhWi + WP-phhy + ShiPh,Ph)- 

This is the desired result (16. 5p and we complete the proof. □ 

Now, we can present a practical postprocessing algorithm which can im- 
prove the accuracy of eigenpair approximations for the Stokes eigenvalue 
problem (12. ip . 

Algorithm 1. 

(1) Solve the discrete Stokes eigenvalue problem /i2.10\) for {uh,Ph,Xh) G 

V;, X Q;, X M. 

(2) Solve the discrete Stokes source problem Ii6. 1\) to get the solution 

{uh,Ph) ^^h^Qh- 

(3) Compute 

The pair {uh,ph,Xh) represent a new (and better than {uh,Ph, Xh)) ap- 
proximation to (u,j9, A). 

Let us discuss two methods to construct the augmented finite element 
space Yh x Qh for solving the Stokes source problem (16.11) . 

Way 1. ("Two grid method" from pBj): In this case, Yh x Qh is the 
same type of finite element space as V/j x Qh on the finer mesh Th with mesh 
size h^{j3 > 1). Here is a finer mesh of Q which can be generated by the 
refinement just as in the multigrid method([35j). 

First, let us consider the case when the exact eigenfunction is smooth and 
have the error estimate (??) and (??). Because the maximum regularity of 
the solution (u,p) of Stokes source problem (I4.19p is {H^{Q))'^ x H'^{Q), we 
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need to chose k < 2. In this case, we obtain the following improved accuracy 
for the eigenpair approximation when (3 = 2([35J) 

\~Xh- X\<Ch^'' forfc<2, (6.7) 

||ii,,-u||i + ||p,,-j9||o<C/i^'' forfc<2. (6.8) 

When f2 is a convex polygonal domain, with the error estimate (??), (??) 
and Theorem 16.11 we have 

\X-~Xh\<Ch^', (6.9) 
\\uh-u\\^ + \\ph-p\\o<Ch^', (6.10) 

where we also choose /? = 2. From the error estimate above, we can find 
that the postprocessing method can obtain the convergence order as same as 
solving the Stokes eigenvalue problem on the finer mesh Th. This improve- 
ment costs solving the Stokes source problem on a finer mesh with mesh 
size 0{h?). This is better than solving the Stokes eigenvalue problem on the 
finer mesh directly, because solving Stokes source problem needs much less 
computation than solving Stokes eigenvalue problem. 

Way 2. ("Two space" method from [30]): In this case, Yh x Qh is defined 
on the same mesh Th but one order higher than V/j x Qh- Since the max- 
imum regularity of the solution (u,p) for the Stokes source problem fl4.19p 
is {H^{VL))'^ X H'^[Q), we can only use the first order finite element space 
to solve the original Stokes eigenvalue problem f l2.10p . and solve the Stokes 
source problem (16. ip in the second order finite element space. So, we only 
have the following error estimate for {uh,Ph, Xh) 

\X-Xh\<Ch\ (6.11) 

\\u-Uh\\i + \\p-ph\\o<Ch, (6.12) 

\\u-Uh\\-i<Ch'^. (6.13) 

First, if the domain Q is smooth, we have the following error estimate 

\X-Xh\<Ch\ (6.14) 

||u - Uh\\i + Hp - Ph\\o < Ch?. (6.15) 

This is an obvious improvement than (16. lip and (16.120 . 

When is a convex polygonal domain, from the regularity of the Stokes 
source problem and the error estimates (16. 3p . (16. 5p and (16.61) . we have 

\X-Xh\<Ch^+^\ (6.16) 
\\vi-^h\\i + \\p-phh<Ch}^^. (6.17) 

This estimate is also an obvious improvement than (16. lip and (I6.12p . 

The improved error estimate above just cost solving the Stokes source 
problem on the same mesh in the second order finite element space. 
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7 Numerical results 



In this section, we give a numerical example to illustrate the efficiency of 
the postprocessing algorithm derived in this paper. Since we do not know 
the exact solution of the Stokes eigenvalue problems, the numerical results 
only give the behaviors of eigenvalue approximations by the postprocessing 
algorithms. 

We consider the Stokes eigenvalue problem (11.11) on the domain Q = 
(0,1) X (0,1). From [34] and ^14j, we choose a sufficiently accurate ffist 
eigenvalue approximation A = 52.3446911 as the ffist true one. 

We ffist give numerical results of the postprocessing algorithm which the 
enriched spaces constructed by refining the current mesh by the regular way. 
Here we use the element (V/i, Dh) = {Pi, P'^^f'^) with 

Pi ={v G H\n) : v\k G Pi{K), \JK G %], 
={v G L\^) : v\k G F-i(K), G T,}, 

to solve the Stokes eigenvalue problem (I2.10p and the Stokes source problem 
(16. ip . The numerical results are shown in Figure [TJ Then we give numerical 



Correction method by refining mesh 

10 I ' 




number of elements 



Figure 1: Errors for refining mesh method with aK = 0.1 

results of the postprocessing algorithm which the enriched spaces constructed 
by one order higher finite element. We first solve the Stokes eigenvalue prob- 
lem (I2.10p by the lowest order stabilization element (V/^, Dh) = (Pi, pfi^c-j 
solve the Stokes source problem (16.11) by second order stabilization element 
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(y,,D,) = (p/,pf-) ([2S]) with 

={v e H\Q) : v\k G P2{K) © • W, VK G T,}, 
P{ii^= ={t; G L2(^]) : ^;|^ e Pi(K), G T,}, 

on the same triangular meshes, where the bubble function tpx is defined 
by the barycenter coordinates Xi^k,^2,k and Xs^k on the element K with 

fK '■= Xi^kX2,kX3^K- 

The numerical results are shown in Figure [2j From Figures [1] and [21 we 



Correction method by higher order finite element 
10^ I • • • 




10"^l — — — — ■ ■ I 

1 ? S 4 S 

10 10 10 10 10 

number of elements 

Figure 2: Errors for higher order method with = 0.1 

can find that the postprocessing algorithm can improve the accuracy of the 
eigenvalue approximations and confirm the theoretical analysis. 

8 Concluding remarks 

In this paper, the LPS method is applied to obtain the approximations of 
Stokes eigenvalue problem and a type of postprocessing method is also pro- 
posed to improve the convergence order for the eigenpair approximation. The 
theoretical analysis is given and the corresponding numerical examples are 
also used to confirm the analysis. The postprocessing method proposed here 
can be coupled with the adaptive mesh refinement in the two-grid method. 
The application of LPS method makes the implementation of adaptive mesh 
refinement more easily for solving Stokes eigenvalue problems especially on 
the meshes with hanging nodes (|32j). 
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In the future, we will extend our postprocessing method to the nonsym- 
metric Stokes eigenvalue problems which is more general in the study of 
linearized stability for the Navier-Stokes equations (|21]). 
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